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ABSTRACT 

The propagator and corresponding path integral for a system of identical par- 
ticles obeying parastatistics are derived. It is found that the statistical weights of 
topological sectors of the path integral for parafermions and parabosons are sim- 
ply related through multiplication by the parity of the permutation of the final 
positions of the particles. Appropriate generalizations of statistics are proposed 
obeying unitarity and factorizability (strong cluster decomposition). The realiza- 
tion of simple maximal occupancy (Gentile) statistics is shown to require ghost 
states. 
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Generalizations of the standard concepts of fermions and bosons have been 
extensively considered in the last few decades. These can be broadly divided into 
"phenomenological" (anyons [1-4], exclusion statistics [5]), which are meant to give 
a better description or understanding of systems of otherwise ordinary fermions or 
bosons, and "fundamental," which would be genuinely new kinds of particles. A 
nice review of these approaches can be found in [6] , and a short and concise resume 
of some relevant results can be found in [7]. 

The first ever consistent extension of fundamental statistics, given by Green 
[8], is parastatistics. In that, the standard bosonic or fermionic fields which would 
create identical particles are replaced by composite fields whose components com- 
mute with themselves and anticommute with each other for parabosons, or vice 
versa for parafermions. The number of components of the fields p defines the "or- 
der" of parastatistics. In general, one can put at most p parafermions in a totally 
symmetric wavefunction, and at most p parabosons in a totally antisymmetric one. 
The degeneracies of occupation of more general multiparticle states are in princi- 
ple calculable but rather complicated. Parastatistics in this approach has been 
well-studied [9-12]. 

The above is a field theoretic realization of parastatistics. Just as in the case 
of fermions or bosons, one can deal with a parastatistical system at a fixed particle 
number in a first-quantized formalism. In this approach, due mainly to Messiah 
and Greenberg [13,14], the iV-body Hilbert space is decomposed into irreducible 
representations (irreps) of the particle permutation group Spj. Since the particles 
are indistinguishable, this group should be viewed as a "gauge" symmetry of the 
system, and states transforming in the same representation have to be identified. 
Moreover, since all physical operators are required to commute with the permuta- 
tion group, each irreducible component is a superselection sector. Therefore, one 
can project the Hilbert space to only some of the irreps of Sjy. Further, only one 
state in each irrep need be kept as a representative of the multiplet of physically 
equivalent states. The resulting reduced space constitutes a consistent quantiza- 
tion of N indistinguishable particles. The choice of included irreps constitutes a 
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choice of quantum statistics. In particular, parabosons correspond to including 
only irreps with up to p rows in their Young tableau, while parafermions to ones 
with up to p columns. Clearly the cases p = 1 reduce to ordinary fermions and 
bosons. 

This description relies on a canonical quantization of the many-body system. 
It is of interest to also have a path-integral formulation of a quantum system, 
since this complements and completes the conceptual framework and usually offers 
orthogonal intuition in several cases. For ordinary statistics this question was 
studied by Laidlaw and DeWitt [15]. In this paper, we provide such a realization for 
parastatistics, or, in general, for any statistics where the Hilbert space is embedded 
in the tensor product of iV one-particle Hilbert spaces (note that this excludes 
anyons and braid statistics). 

The starting point will be the coordinate representation of the full (unpro- 
jected) Hilbert space, spanned by the position eigenstates \x\,...xn >= \x > 
(where be in a space of any dimension). The collection of such states for a 

set of distinct X{ transforms in the iV!-dimensional defining representation of Sjy 

P\x >= \Px >= |xp-i(i), • • • Xp-i(pf) > (1) 

where P is a permutation (the appearance of P _1 in the above is necessary so 
that products of permutations be represented in the right order). If any of the 
coordinates x% coincide the above is not the full defining representation any more. 
The set of such states, however, is of measure zero (the coordinate space assumed 
continuous) and thus they can be safely ignored. (We assume that there are no 
interactions singular at coincidence points that might dynamically make such states 
of nonzero weight.) 

Projecting the Hilbert space to an irrep R of Sn amounts to keeping only linear 
combinations of states within this multiplet transforming in R, that is, 

\a;x >= y ^2 / C a (P)P\x > , a = l,...d R , d R = dim(R). (2) 
P 
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where the sum is over all elements of the permutation group and C a (P) are appro- 
priately chosen coefficients. If we denote with R a (,(P) the matrix elements of the 
permutation P in the representations R, 

P\a,x>=J2 R ab(P~ 1 )\b,x> (3) 

b 

The defining representation decomposes into irreducible components, classified by 
Young tableaux, each appearing with a certain multiplicity. Should we keep only 
one irrep out of each multiplicity or the whole multiplet? To decide it, note that 
if instead of the base state \x > for the construction of the states \a,x > we 
choose a different permutation P \x >, then although the new states \a,P x > 
constructed through (2) still transform in the irrep R, in general they are not 
linear combinations of \a, x > but rather span a different copy of R. Since we can 
continuously move in the configuration space from \x > to P Q \x >, we conclude 
that we must keep a//irreps R within each multiplet. (In other words, although for 
each point in the Hilbert space \x > this multiplet is reducible, the fiber of these 
representations over the Hilbert space is connected and irreducible.) 

To realize explicitly the above, we construct the states 



ab, x > = 




Using the group property of the representation R(Pi)R(P2) = R(P\P2), we deduce 
that under the action of the group Sn and under change of base point x the above 
states transform as: 

P\ab,x >= ^R ac (p- l )\cb,x > , \ab, Px >= R cb {P~ l )\ac, x > (5) 

c c 

Thus we see that the first index in these states labels the different elements of a 
single irrep R, while the second index labels the different equivalent irreps in the 
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multiplet. Since both indices take cIr values, we recover the standard result that 
each irrep of Sn is embedded in the defining representation a number of times 
equal to its dimension. 

Consider now the matrix element < ab, x\A\cd, y >, where A is any physical 
operator, that is, any operator commuting with all elements P of SV- Substituting 
the definition (4) and using the unitarity of P (pt = p-i) an d of R (R* ab (P) = 
Pfea(P _1 )) we obtain, after a change in summation variable, 

< ab, X \A\cd, y>=^ J2 Rbe(P')Rea(P- l )Rcd(P) < x\AP'\y > (6) 

' P,P',e 

Using further the orthogonality (Shur's) relation (see, e.g., [16]) 

V P ab (P)P crf (p- 1 ) = — 5 a Ac (7) 
p 

we finally obtain 

< ab, x\A\cd, y>=J2 S acRbd(P) < x\A\Py > (8) 
P 

Let us first choose A = 1. Then the above provides the overlap between the states 

< ab, x\cd, y>=Y^ S ac Rbd(P)S(x - Py) (9) 
P 

For x in the neighborhood of y it is P = 1 which contributes to the normalization, 
for which RbdiX) — $bd an d we recover the standard continuous normalization 
between the states. 

Now we can choose A = e~ , where H is the hamiltonian, and thus find the 
propagator G(ab, x; cd, y\t) between the states of the system. It is clear from (8) 
that the first index a in the state \ab,x > propagates trivially. Since this is the index 
that corresponds to the different but physically equivalent states within each irrep 
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R, we conclude that the required projection of the Hilbert space to the physical 
subspace amounts to simply omitting this index from all states. (That is, freeze 
this index to the same fixed value for all states of the theory; no physical quantity 
will ever depend on the choice of this value.) On the other hand, the second 
index, corresponding to different equivalent irreps, does not propagate trivially 
and must, as argued before, be kept. We are led therefore to the physical states 
1 6a, x >— > \a,x > and the propagator 

G R (a, x; b,y\t) = ^ R ab (P) G(x, Py; t) (10) 
P 

where G{x, Py; t) = < x\e~ iHt P\y > is the usual many-body propagator. We note 
that, due to the transformation property (5), the states \a,Px > are linear com- 
binations of states \a,x >. Therefore, projecting down to the physical subspace 
corresponding to R amounts to trading the original N\ copies of physically equiv- 
alent states \Px > for a number dp of global internal degrees of freedom for the 
system, labeled by the index a. 

It is now easy to write down the path integral corresponding to identical par- 
ticles quantized in the i?-irrep of Sjy. G(x, Py; t) can be expressed as an iV-body 
path integral in the standard way, with particles starting from positions X{ and 
ending in positions Pyi = yp-i^y Since all permutations of particle positions are 
physically equivalent, (10) instructs us to sum over all sectors where particles end 
up in such permuted positions, weighted with the factors R a t,(P) depending on the 
internal degrees of freedom of the initial and final states. From (4), (9) we can 
write the completeness relation within the physical subspace 

Ir= f ^—f^2\a,x >< a,x\ (11) 

J ' a 

and with the use of (11) it is easy to prove that the above path integral is unitary, 
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that is, 

J ' b 

The extension to parabosons, parafermions or any similar statistics is immedi- 
ate. Let S = {Ri, . . . R n } be the set of allowed irreps of Sn in the Hilbert space. 
The internal degree of freedom now takes values A = (R,a), where R E S and 
a — 1, . . . dji labels the internal degrees of freedom within each irrep. So, overall, 
A takes + ■ ■ ■ dR n different values. The propagator (and corresponding path 
integral) is obviously 

G s (A,x;B,y\t) = ^ S(P) AB G(x, Py;t) , where S(P) AB = 5 RA ,R B (R A ) ab (P) 
P 

(13) 

For parabosons (parafermions) of order p, S is the set of Young tableaux with up 
to p rows (columns). We note that the irreps for parafermions are the duals of 
those for parabosons (the dual of a tableau is the tableau with rows and columns 
interchanged). In an appropriate basis, the representation matrices of dual irreps 
R, R are real and satisfy 

R ab (P) = (-l) p R ab (P) (14) 

where (— l) p is the parity of the permutation. We arrive then at the relation 
between the weights for parabosons and parafermions of order p: 

S p f(P)ab = (-1) P S pB (P)ab (15) 

This extends a similar relation for ordinary fermions and bosons, for which there 
are no internal degrees of freedom and Sb(P) = 1- 

From the path integral we can evaluate the partition function, by simply shift- 
ing to the euclidean periodic propagator Ge{(3) = e~@ H and summing over all 
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initial and final states, with the measure implied by (11). Given that 



J2Raa(P)=tTR(P)=XR(P) 



(16) 



a 



we get the expression in terms of the characters of Sn 



Zs{T)= —Y. S ( P ) <x \ G E^)\ Px >i where S(P) = Y,XR(P) (17) 



The interpretation in terms of a periodic euclidean path integral is obvious. The 
characters xr(P) are a se t °f integers, and thus the "statistical factors" S(P) 
weighing each topological sector of the path integral are (positive or negative) 
integers. In the case of parabosons of any order p, however, we note that the 
statistical weights are positive (or zero) integers. The ones for parafermions can be 
either positive or negative, as given by 



We do not have a general formula for S p b(P) for arbitrary p. 

From the above results we can derive the partition function for a gas of parasta- 
tistical particles as well as the allowed occupancy of single-particle states. Consider 
a collection of non-interacting particles, for which the hamiltonian is separable into 
a sum of one-body hamiltonians H = J^i^^i)- Let the energy eigenvalues of 
the one-body problem be ej and the corresponding one-body Boltzmann factors 
Zi = e~ /3ei . Consider now a sector of the euclidean path integral characterized by 
the permutation of final points P. It is clear that this path integral Zp decomposes 
into a product of disconnected components, characterized by the fact that the par- 
ticle worldlines in each component mix particles only within the same component. 
This means that, within each component, particles mix under a cyclic permutation 
(since following the worldline of each particle must successively lead to every other 



P 



Res 



S pF (P) = (-l) p S pB (P) , S pB (P) > 



(18) 
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particle in the component). Each element P of Sjq is then decomposed into a prod- 
uct of commuting cyclic permutations. The number of particles n participating in 
each cyclic permutation constitute the cycles of P (obviously Yl n = N). We are 
led, thus, to the fact that for noninteracting particles 

z P = n z n (19) 

necycles(P) 

The path integral Z n for a cyclic permutation of n particles, on the other hand, 
can be thought of as the path integral of a single particle winding n times around 
euclidean time (3. This means that 

Z n {(3) = Zi (n/3) = J2 z i= W " N ( 2 °) 

i 

and the corresponding expression for Zp becomes 

Z P = J] W n [ Zi ] (21) 

necycles(P) 

The expression for the full partition function then becomes 

z «=eem^) n w *M ( 22 ) 

ReS P ' necydes(P) 

We recognize the sum over P in (22) as Frobenius' relation, connecting the sum 
over the Schur functions W n [zi] to the characters of SU(M) xr\ z i\- We get the 
final result 

Zs=E ^l = E d ^ ) (23) 

ReS ReS ael \ z i I 

where £j is the length of the j-th row of the Young tableau of R. This reproduces 
the result of Chaturvedi for the partition function [17] and an earlier result of 
Suranyi for p = 2 [18]. We stress that the above result holds only for noninteracting 



9 



particles. For interacting particles the topologically disconnected components Z n 
of the path integral Zp are still dynamically connected and factorization fails. One 
has to go back to the full expression (17) for the partition function in that case. 

To find the degeneracy of states, we need to decompose X-rN] appearing in 
(23) in monomials Yii Z T 

£x*N = I>MlK° ( 24 ) 

R&S { Pl } i 

The coefficients D[pi] of these monomials, called Kostka-Foulkes numbers [17, 19], 
are non-negative integers which determine the degeneracy of the state with pi 
particles occupying each energy level €{. To find these integers in a systematic 
way, we use the following trick: consider that the particles are bosons and have 
an internal degree of freedom transforming in the fundamental of SU(N) (in fact, 
SU(M) with M > N would also do). Since under total permutation of particle 
coordinates and internal degrees of freedom the states must transform trivially, we 
conclude that the irrep of the color SU(N) for each state must be the same as the 
irrep of the coordinate permutation group Sn (meaning they have the same Young 
tableau and thus the same symmetries). A state with p particles in the same 
level €i transforms in the p-fold symmetric irrep of SU(N). Therefore, a state 
with occupancies pi for each level transforms under the direct product of pj-fold 
symmetric irreps, one for each level e\. Decomposing this product into irreducible 
components, we will obtain each representation R of SU(N) a number of times 
Dji[pi\. Each such irrep will transform under a similar irrep of Sn and thus will 
correspond to a unique physical state in the quantization of N identical particles in 
the representation R. Therefore, the degeneracy D[pi] can be found by summing 
the number of times that each allowed irrep R G S appears in the direct product 
of symmetric irreps pi, which can be found using standard SU (N) Young tableaux 
composition rules. We also see that, if the internal degree of freedom group is 
chosen to be SU (p) (where p may be smaller that N), we will only get irreps with 
up to p rows. Therefore, we recover the known result that parabosons of order p 
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can be viewed as bosons with an internal SU (p) symmetry, where we identify each 
irrep of SU(p) as a unique physical state [11,12]. A similar construction can be 
repeated starting from fermions instead of bosons. We recover a dual expression 
for the degeneracies D[pi], where now we form the direct product of p^-fold fully 
antisymmetric irreps of SU(N), and a similar expression of parafermions of order 
p as fermions with an SU (p) internal symmetry. 

As was argued in [10,14], parastatistics particles obey the cluster decomposition 
principle, in the sense that the density matrix obtained by tracing over a subset of 
particles which decouple from the system can be constructed as a possible density 
matrix of the reduced system of remaining particles. From (17), however, we see 
that the partition function of two dynamically isolated sets of particles N\ and 
N2 does not factorize into the product of the two partition functions, since the 
statistical weights S(P) in general do not factorize into S(Pi)S(P2) when P is 
the product of two commuting elements Pi and P2. Equivalently, this means that 
the occupation degeneracy D[pi] does not factorize into the product of individual 
occupation degeneracies for each level ej. This has important physical implications. 
If the two sets of particles are totally isolated, it does not make sense to evaluate the 
partition function of the total system, since the statistical distribution can never 
relax to the one predicted by that partition function. The individual partition 
functions of the subsystems are the relevant ones. If, however, the two sets are 
only weakly coupled, then initially each set will distribute according to its reduced 
partition function, but after some relaxation time (depending on the strength of 
the coupling between the two sets) they will relax to the joint distribution function, 
which, we stress, will not even approximately equal the product of the individual 
ones. Thus, cluster decomposition holds in an absolute sense but fails in a more 
realistic sense. In contrast, fermions and bosons respect cluster decomposition in 
both senses. 

The obvious generalization of quantum statistics, based always on the assump- 
tion that the many-body Hilbert state is embedded into the tensor product of 
many one-body Hilbert spaces, is to generalize the set of allowed irreps S beyond 
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the one relevant to parastatistics. We may, however, further include more than one 
state for each included irrep of S]y. This seems unmotivated, in view of the fact 
that such states are physically indistinct, but it is certainly consistent. It could 
mean, for instance, that the particles have some hidden internal degrees of freedom 
accounting for the extra degeneracy, which are invisible to the present hamiltonian 
but may become dynamically relevant later. The most general situation, then, is 
that we include Cr states from each irrep R. The generalization of all previous 
formulae for this case is quite immediate, S(P) and Cr being related by 

s{p) = Y,Crxr{p), Cr = ^J2 s{p)xr{p) (25) 

R ' P 

The case of distinguishable particles ("infinite statistics" [12,20,21]), in particular, 
is reproduced by accepting all states in each irrep, that is, Cr = d,R. Since R 
appears exactly d,R times in the defining representation of Sn, S(P) above be- 
comes the trace of P in that representation. But all P ^ 1 are off-diagonal in 
the defining representation, so we get Si n f(P) = N\8p t i, recovering the standard 
distinguishable particles result. 

We summarize here by pointing out that the most general statistics of the type 
examined here is parametrized by any of three possible sets of numbers. The first 
is, as just stated, the number of states Cr accepted for each irrep R of Sn. Since 
the irreps of Sn are parametrized by the partitions of N (lengths of rows or the 
Young tableau), there are as many Cr as there are partitions of N. The second set 
is the statistical weights S(P) appearing in the partition function (euclidean path 
integral). Clearly these weights are invariant under conjugation of P — > QPQ^ 1 , 
since this simply amounts to a relabeling of the particle worldlines. Thus S(P) 
depends only on the conjugacy class of P, that is, the cycles of P. The possible sets 
of cycles are the same as the partitions of N; so, again, the S(P) are numbered by 
partitions of N. Finally, we could use the degeneracy of a many-body occupancy 
state D[pi] as our definition. There are as many ways to distribute particles in 
one-body states as there are partitions of N, so this set also has the same number 
of elements as the previous two. 
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What are the restrictions or criteria to be imposed on the above parameters? 
The first one is unitarity, that is, the existence of a well-defined Hilbert space with 
positive metric. This requires that Cr be non-negative (no negative norm states) 
integers (no "fractional dimension" states). The other will be what we call "strong 
cluster decomposition principle," that the partition function of isolated systems 
factorize. This is a physical criterion, rather than a consistency requirement. To 
summarize: 

• Unitarity : Cr non — negative integers 

• Strong cluster decomposition : S(P) = 1 | S(n) or D[pi] = TT D(pi) 



The strong cluster decomposition, in particular, implies the existence of a grand 
partition function, obtained (in the case of noninteracting particles) by exponen- 
tiating the sum of all connected path integrals (P a cyclic permutation of degree 
n) with weights S{n)/n (1/n is the symmetry factor of this path integral, corre- 
sponding to cyclic relabelings of the particles). The grand partition function will 
further factorize into a product of partition functions for each level €{. Thus, S(n) 
are cluster coefficients connected to D{n) in the standard way 



The above formula, in fact, provides the easiest way to relate D[pi] and S(P) in 
the general case (no strong cluster property): simply expand the right-hand side 
of (26) in powers of z and substitute every term S(ni) ■ ■ ■ S(nk) with S(ni, . . . nk). 
This gives D(p). To find D[pi] = D(pi, . . .p^) simply evaluate D(p\) ■ ■ ■ D(pk) 
using the above formula and again consolidate each product S{n\) • • • Sin^) into a 
single S(ni,...n k ). 

If we assume that S(l) = D(l) = 1, then it is easy to verify that the only 
solution of the above two criteria is ordinary fermions and bosons. The situation 
is different, however, when S(l) = D(l) = q > 1 (this would mean, e.g., that the 



n£cycles(P) 




(26) 
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particles come a priori in q different "flavors"). The possibilities are manifold. All 
these generalized statistics share the following generic features: 

• The degeneracy of the state where n particles occupy different levels is q n . 
(Indeed, D(l, 1, ... 1) = D(l) n = q n .) 

• If state A can be obtained from state B by 'lumping' together particles that 
previously occupied different levels, then D(A) < D(B). (E.g., D(3) < D(2, 1) < 

£>(i,i,i).) 

This second property is actually related to the (weak) cluster decomposition as 
formulated in [14], which is obviously covered by the strong cluster property. 

The above possibilities include the obvious special cases of q\ bosonic flavors 
and q2 fermionic ones (q\ + qi — q), for which S(n) — q\ — {—l) n q2 ) along with 
many other. As an example, we give the first few degeneracies for many-particle 
level occupation for all statistics with q = 2: 



£>(1) 


= 2, D{2) 


= 4, D(3) 


= 8 


D(l) 


= 2,D(2) 


= 3, D(3) 


= 6,5,4(5 + 5) 


D(l) 


= 2,D(2) 


= 2, D(3) 


= 4,3,2(5 + F),l,0 


D(l) 


= 2,D(2) 


= 1,D(3) 


= 2,1,0(F + F) 


D(l) 


= 2,D(2) 


= 0,D(3) 


= 



The specific choices denoted by B + B, B + F and F + F are the ones cor- 
responding to two bosonic, one bosonic and one fermionic, and two fermionic fla- 
vors respectively. The topmost statistics could be termed "superbosons" and the 
bottom one "superfermions" of order 2. We also remark here that the "(p, q)- 
statistics" introduced in [14] can be realized as particles with p bosonic and q 
fermionic flavors, where we identify each multiplet transforming irreducibly under 
the supergroup SU (p, q) as a unique physical state. 

Finally, we direct our attention to the first known attempt to generalize the 
ordinary Fermi or Bose statistical mechanics, by Gentile [22]. The rule is simply 
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that up to p particles can be put in each single-particle level. This corresponds 
to D{n) = 1 for n < p, and D(n) = otherwise. This has been criticized [7] 
on the grounds that fixing the allowed occupations for each single-particle state 
is not a statement invariant under change of single-particle basis. It is clear that, 
in the language of this paper, any statistics satisfying the unitarity requirement is 
consistent and basis-independent. Therefore, Gentile statistics must violate uni- 
tarity. Indeed, it is easy to check that all weights Cr for such statistics are integers 
(this is generic for all statistics with integer D(n)), but not necessarily positive. 
In the specific case of p — 2, e.g., where up to double occupancy of each level is 
allowed, the degeneracies of each irrep of Sn (parametrized, as usual, by the length 
of Young tableau rows) up to iV = 5 are 

C 2 = C 2 i = C22 = C221 = 1, Cm = Cini = C2111 = -1, else C R = (27) 

We see that representations 111, 1111, 2111 correspond to ghost (negative norm) 
states and their effect is to subtract (rather than add) degrees of freedom. We 
also remark that the path integral realization of exclusion statistics exhibits both 
negative and fractional statistical weights, signaling breakdown of unitarity [23]. 
This is inconsequential in that case, since exclusion statistics is valid only as a 
macroscopic (statistical) description of some (interacting) systems of particles. 

In conclusion, we have presented the many-body propagators and correspond- 
ing path integrals of particles obeying parastatistics or any other type of statistics 
based on irreps of the permutation group. We argued that there are many pos- 
sible unitary generalizations obeying the strong cluster decomposition principle, 
although they all require more than one flavor of particles. Several other direc- 
tions of investigation and open questions suggest themselves. To name a few, the 
statistical mechanics of such generalized statistics particles should be examined. 
Also, it should be checked if they can be realized as particles with specific hid- 
den internal symmetries and an appropriate projection of the Hilbert space, in a 
fashion similar to parastatistics. Independently, it would be interesting to see if 
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Gentile statistics can be consistently realized by introducing 'benign' ghosts which 
account for the negative norm states while decoupling from all physical processes, 
just as in gauge theories. Finally, a similar analysis could be attempted for gener- 
alized particle statistics in 2 + 1 dimensions. In fact, a similar treatment, based on 
the permutation group, has been used to obtain perturbative results for anyonic 
particles [24-26] . It would be interesting to examine whether non-abelian irreps of 
the braid group, instead of the permutation group, could be considered. 

Acknowledgements: I would like to thank J. Myrheim for discussions. 

REFERENCES 

1. J.M. Leinaas and J. Myrheim, Nuovo Cimento 37B, 1 (1977). 

2. G.A. Goldin, R. Menikoff and D.H. Sharp, J. Math. Phys. 21, 650 (1980); 
22, 1664 (1981); Phys. Rev. D28, 830 (1983). 

3. F. Wilczek, Phys. Rev. Lett. 48, 1114 (1982); 49, 957 (1982). 

4. Y.S. Wu, Phys. Rev. Lett. 52, 2103 (1984). 

5. F.D.M. Haldane, Phys. Rev. Lett. 67 (1991) 937. 

6. J. Myrheim, Anyons (Notes for the Course on Geometric Phases, ICTP, 
Trieste 6-17 Sept. 1993). 

7. O.W. Greenberg et al, [hep-ph/ 9306225 . 



8. H.S. Green, Phys. Rev. 90, 270 (1953). 

9. O. Steinmann, Nuovo Cimento 44, A755 (1966). 

10. P.V. Landshoff and H.P Stapp, Ann. of Phys. 45, 72 (1967). 

11. Y. Ohnuki and S. Kamefuchi, Phys. Rev. 170, 1279 (1968); Ann. of Phys. 
51, 337 (1969). 

12. S. Doplicher, R. Haag and J. Roberts, Comm. Math. Phys. 23, 199 (1971); 
35, 49 (1974). 



1G 



13. A.M.L. Messiah and O.W. Greenberg, Phys. Rev. B136, 248 (1964); B138, 
1155 (1965). 

14. J.B. Hartle and J.R. Taylor, Phys. Rev. 178, 2043 (1969); R.H. Stolt 
and J.R. Taylor, Phys. Rev. Dl, 2226 (1970); Nucl. Phys. B19, 1 (1970); 
J.B. Hartle, R.H. Stolt and J.C. Taylor, Phys. Rev. D2, 1759 (1970). 

15. M.G.G. Laidlaw and CM. DeWitt, Phys. Rev. D3, 1375 (1971). 

16. M. Hammermesh, Group Theory, Addison- Wesley Eds. (1962). 

17. S. Chaturvedi, U. Hyderabad preprint, hepth/9509150. 

18. P. Suranyi, Phys. Rev. Lett. 65, 2329 (1990). 

19. C. Kostka, Crelles Journal 93, 98 (1882); H.O. Foulkes, Permutations, 
Eds. Gauthier-Villars, Paris (1974). 

20. A.B. Govorkov, Theor. Math. Phys. 54, 234 (1983). 

21. O.W. Greenberg, Phys. Rev. Lett. 64, 705 (1990). 

22. G. Gentile, Nuovo Cimento 17, 493 (1940). 

23. A.P Polychronakos, Phys. Lett. B365, 202 (1996). 

24. J. Myrheim and K. Olaussen, Phys. Lett. B299, 267 (1993). 

25. J. Desbois, C. Heinemann and S. Ouvry, Phys. Rev. D51, 942 (1995). 

26. A.D. de Veigy, Nucl. Phys. B458 [FS], 533, 1996). 



17 



